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Abstract 



We apply the retarded/advanced formalism of real time field theory to the QED or 
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I. Introduction 

In the study of field theory at high temperature [1] , the development of resummation 
techniques [2-4] has considerably extended the domain of validity of the perturbative ap- 
proach. From the practical point of view, several possible signals for the formation of the 
quark-gluon plasma have now been studied in this framework. In particular the production 
of soft virtual photons [5,6] as well as hard real photons [7,8] have been discussed in great 
detail. Besides providing a consistent perturbative expansion, the resummed series gives 
a finite result for quantities which would otherwise not be well defined in a more naive 
approach. 

The resummation technique of "hard loop" expansion has been formulated in the 
imaginary-time formalism (ITF) [1,9,10] of thermal field theories. On the other hand, 
there has been many developments in the real-time formalism (RTF) [1,10-12] or in the 
thermo-field dynamics (TFD) approach [13]. The latter seems more appropriate for the 
study of dynamical phenomena since time is kept as an independent variable throughout 
the calculation. Equivalently, in momentum space, no analytic continuation has to be per- 
formed to obtain physical quantities. Several studies have been devoted to the comparison 
of the ITF and the RTF [14-20]. In particular, a simple way to relate the two formalisms 
consists in constructing in the RTF the retarded/advanced (R/A) Green's functions which, 
at least for the two- and three-point functions, have a close connection to the ITF Green's 
functions. In this paper, we pursue the study of the R/A functions and show, in QED, 
how to formulate the hard-loop expansion: several formulae are derived which are useful 
in the calculation of physical quantities. 

In the R/A formulation of the RTF, the propagators are still 2x2 matrices but they 
are defined to be diagonal and constructed, at least at the lowest order of perturbation 
theory, from the retarded and advanced propagators of the T = theory [17,18]. As 
for the vertices, they become temperature dependent. The R/A (amputated) Green's 
functions are in fact related to the causal combinations of RTF Green's functions first 
introduced by Kobes [14,15]. In TFD, which has very similar Feynman rules to the RTF, 
such a diagonalization can also be performed and it is related to a thermal Bogoliubov 
transformation, thereby allowing the interpretation of the diagonal propagators as those of 
statistical quasi-particles [21]. This formulation seems therefore quite naturally appropriate 
to the study of thermal field theories. 

In the following, we first complete the study started in [17] by constructing explicitly 
the diagonalization for fermion propagators and defining the Feynman rules in the R/A 
approach. We then discuss, in QED, the one loop expression for the 2-, 3- and 4-point 
functions in the context of the hard loop approximation and we easily recover the ITF 
results continued to real energies. We then set up a general formula, at the multi-loop level, 
for the 2-point functions and discuss how it is related to the ITF results. As an application, 
we consider the production of soft real photons and show that, despite cancellations due to 
thermal gauge invariance between various terms of the hard loop expansion, there survives 
a divergence. One of the advantages of the R/A approach over the ITF one is that it does 
not require any analytical continuation. Furthermore the spin structure of the diagrams is 



1 



the same as at T = and no continuation from Euclidian to Minkowski space is needed. 
Also, the possibility of using contour integration may simplify the calculation. 



II. The propagators and vertices in the R/A formalism. 

In the real-time formalism, the fermion propagator, defined on a contour characterized 
by a (fig. 1) can be written as a product of 2 x 2 matrices [17]. 

S F (P) = (/> + M) U F (P) D(P) V F (P). (1) 

Denoting the retarded and advanced bosonic propagators at temperature 

1 1 



P 2 — M 2 + iepo 2Q \p - Q + ie p + Q + ie 
i i ( 1 1 



(2) 



P 2 — M 2 — iepo 20 \po - Q-ie p + Q-ie 
(O = (p 2 + M 2 ) 1 / 2 ) we introduce the diagonal matrix 

D( P) =(zu = ( A f> A ; (P) ) (3) 

The other matrices will be specified shortly after we define the photon propagator. In the 
Feynman gauge, the latter takes the form 

G^(P) = -g^ U B (P) D(P) V B (P) (4) 

where the same diagonal matrix as in eq. (1) appears. The diagonalization matrix, denoted 
generically are defined in a similar fashion for bosons and fermions, and besides 

the contour a, they depend on arbitrary scalar functions b(p), c(p) 

UM(P) = (t© = -„ nM(- P0 ) ( 6 _?;_ I „ po V*- 1 *-** ) ( 5 ) 

where, not surprisingly, rj = 1 ([rj\ = B) for a boson and rj = — 1 ([rj\ = F) for a fermion 
and n^(po) is the usual Bose-Einstein or Fermi-Dirac distribution. f3 is the inverse of the 
temperature. The chemical potential has been set to 0, but it is easily introduced by the 
substitution (3po —> (3(po — p), crpo remaining invariant. 

The basic principle, in the R/A formalism, is to associate the diagonalization matrices 
to the vertices in a natural way, keeping therefore the matrices D as diagonal propagators. 
This leads to different types of vertices depending on the momentum flow. For instance, 
for all incoming momenta as in fig. 2a (P + Q + R = 0) we introduce 

-H a ps(P Q, R) = -ig abd V F a (P) V*(Q) V[ d {R) (7) 
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The latin indices refer to the 1 (particle) or 2 (ghost) fields of the usual formulation of 
RTF so that gm = e, (7222 = — e (e is the electron charge), all the other couplings being 
0. The greek indices take the values R or A. In defining the vertex function 7,^,5 we leave 
out the Dirac or Lorentz part of the vertex and keep only its scalar R/A structure. From 
the definitions above it can be shown that 

laP5 (P, Q,R) = e (b(P)) 5 ^ (b(Q)) s ^ (b(R)) 5 ** 

(-C(P))^ (_ C (Q))<5/3A (- C (R))6sA e aL ^ _ (_ 1 )<5 Q K+5^ e -/3L ) 

with L Q = p 5 aR + q 5f3R + r Q 5 SR . 

It is immediately clear at this point that ^aaa always vanishes while ^rrr = 
because of momentum conservation. These results reflect the causality requirement that 
three particles propagating forward in time (or backward in time) cannot annihilate into 
(or be created from) the vacuum. It is convenient, and perhaps more natural to introduce 
the vertex in figure 2b, where the flow of momentum follows the fermion line. Momentum 
conservation now reads P + Q = R. We define 

-il*pAP, Q> R) = -Wabd V a F a (P) V*{Q) U F (R) (9) 
An expression, similar to eq. (8), but not so symmetrical can be derived 

la p, d {P,Q,R) = -e (b(P)) s "« {b{Q)) hR (b(R))- SsR 

(-c{P)e- apo ) 5 <* A {-c{Q)e- aqo) ) 5 ^ A {-c{R)e- ara )- &SA ( 10 ) 
n F (ro)e^ r ° S5R (— l) SsR + (_iyaR, e -P(po5 a R+qo5pR-r 5 SA )~^ 

Admittedly, this expression is not very illuminating, but it immediately allows, by com- 
parison with eq. (8) to derive the relations 

7 a/9;fl(P, Q; R) = - h ; m ° m i«pa{p, Q, -R) 

b(K)c(-H) 

n F (rn)e° r ° 

7a /3 ;A(P, Q; R) = - b{ _^ c(R) 7a/9fl(P, Q, ~R) 

The obvious choice 

b(-R) c(R) = -n F (r ) e ar ° (12) 

simplifies the "crossing" relation when the fermion momentum R is changed to —R, since 
we can then simply write 

7 a/3;5 (P, Q; R) = >y af) s(P, Q, -R) (13) 

where we have introduced 5 = A, R the conjugate index of 5 = R, A. Similarly, if the 
"crossing" property of the photon line is studied, we find expression such as eq. (11) with 
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the difference that n B (— q ) appears rather than — n F (— q ). The choice, in the bosonic 
diagonalization matrices 

b(-R) c(R) =n B (q ) e aqo (14) 

allows to keep the same crossing property eq. (13) for both fermion or boson lines. Further 
choices can be made to simplify the expression of the thermal vertices. For example, taking 
b = 1, we immediately obtain 

1rra{P, Q, R) = 1arr{P, Q, R) = 1rar{P, Q,R) = e 

lRAA (P, Q, R) = -e ^ (r .° } = -e (1 + n B (q ) - n F (r )) 

(TDn n\ n F (q )n F (r ) . F( , F( ^ 

7ARA\P,Q,R) = -e — W( ■ — = -e 1-n (q ) - n r )). 

(qo + ro) 

Together with 

1RRR = 1 AAA =0 (16) 

and the temperature independent propagators 



MP)^+M)(*f> A ; (p) ) 



(17) 



the diagrammatic rules for thermal QED in the R/A formalism are completely specified. 
Let us note that all reference to the arbitrary contour a has now disappeared from the 
Feynman rules. Unlike eqs. (16), the results eqs. (15) and (17) are not universal. For 
instance, injl8] the propagators are chosen to be anti-diagonal so that reversal of momen- 
tum yields D(—P) = (D(P)) T . Here we keep the diagonal form as it appears more natural 
for the construction of loop diagrams. Furthermore, by a different choice of the arbitrary 
functions b and c, we could have imposed that the vertices of type ^raa be equal to the 
coupling e. The vertices such as ^rra would then have become linear in the statistical 
weights. 

Any (amputated) n-point Green's function in the R/A formalism can be constructed 
from the corresponding ones in the RTF by a generalization of eqs. (7) and (10). The R/A 
functions then appear as specific linear combinations of the usual real-time functions. It 
is more practical, however, if one is interested in higher order perturbative calculations to 
construct the R/A functions directly by application of the rules given above. Despite a 
rather cumbersome notation, the Green's functions are nicely expressed in terms of tree 
diagrams. 

The diagonalization approach to RTF is more than an algebraic trick. The R/A func- 
tions are "more natural" than the Green's functions in the real-time formalism. For exam- 
ple, it has been shown that the imaginary part of the 2-point function (there exists only two 
such functions related by complex conjugation) is proportional to the opacity factor in the 
Boltzman equation for a distribution near equilibrium [14,18]. Furthermore, the n-point 
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functions with all but one indices set equal to R are the causal Green's functions introduced 
by Kobes [14,15]. A different perspective comes from thermo-field dynamics. The Feynman 
rules in TFD and RTF are very similar (for systems in equilibrium which is the case consid- 
ered here). Recently Henning and Umezawa [21] studied the diagonalization of the 2-point 
function in TFD. The diagonalization matrices (Bogoliubov transformations) have two pa- 
rameters: a which characterizes the thermal vacuum and s, with the relations to our free 
parameters given by a = l — a/(3 and s = Inb^l + n(po) = — uica/1 + n(po). In TFD, the 
Bogoliubov matrices parametrize the transformations which take the physical (point-like) 
particles to the (thermal) quasi-particle states and the diagonal (un-amputated) 2-point 
function D(P) describes the propagation of the quasiparticles in the medium. We note also 
that by a differerent Bogoliubov transformation we could introduce the diagonal matrix 
constructed with the usual Feynman propagator and its complex conjugate. 

We turn now to the study of loop diagrams. 



III. R/A functions in the one-loop approximation 



Following the method introduced for \<p 3 we construct 2-, 3- and 4-point functions 
perturbatively. We implicitly assume that we deal with QED or QCD. In general, we do 
not specify the spin structure of the diagrams since it factorizes from the R/A functions 
and it is the same as at temperature. We first consider a generic two-point function with 
external momentum Q and index (3 flowing through the diagram. The internal momenta 
are P, with fermion of boson statistics r)p, and R, with statistics tjr, such that P + Q = R 
[fig. 3]. Applying the rules eqs. (15)-(17) and using the crossing relation eq. (13) we arrive 
at 



/d n P 
— p >y afa (P,Q,-R) A 5 (R) ld(0s (-P,-Q,R) A a (P) 

(l + Vpn [ ^ ] (po)) (A fl (P) - A A (P)) A^R) (lg) 
+ d+ W-Vo)) (A R (R) - A A (R)) Ap{P) 



. / d n P . 



The symbol p denotes the Dirac and spin structure of the diagram. For example, for the 
fermion self-energy we have, in Feynman gauge p = —7^ (/? + M) 7^ if R is the fermion 
internal line. In arriving at the final form of eq. (18) we have dropped in the integrand 
terms of type Ar(P)Ar(R) or Aa(P)A a (R) independent of the statistical weight. These 
terms have poles in the po complex energy plane only on one side of the real axis. By 
closing the po integration contour in the other half-plane they are seen to give a vanishing 
contribution to the 2-point function. By the same token, the factor ^ can be replaced 
by any constant because shifting the numerical value only generates irrelevant A^A^ or 
AaAa factors. We prefer keeping the | factor since, using the relation 

A R (P) - A A (P) = 2n e(p ) S(P 2 - M 2 ), (19) 
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we find that it leads to an invariant form under the reversal of the sign of momentum P: 



Q + ripn^iP)) (A B (P) - A A (P)) = 2tt Q + i/pn^dpol)) <5(P 2 - M 2 ). (20) 

When evaluating eq. (18) we can either use the 5-function constraint or close the Po 
contour in a conveniently chosen half-plane. For instance for the first term in r^(Q) 
encircling poles in the upper half-plane it is seen that only the poles of Aa(P) contribute. 
In particular, the poles on the imaginary axis associated to the statistical function never 
contribute since they have a vanishing residue as A^(P) — Aa(P) = on the imaginary 
axis. For the second term, it is preferable to close the contour in the lower-half plane to 
retain only the singularities of Ar(R) = Ar(P + Q). 

The two-point function obeys some general relations [17,18]. We have 

K a (Q) = T^(Q) (21) 

provided, it is assumed that the complex conjugation does not operate on the spinor 
structure p but only on the R/A part of the integrand. It is based on the property 

A*(P) = -A a (P). (22) 

We can also prove, for QED/QCD like theories and assuming massless particles, that 

r aa (-Q) = ±r aa (Q) (23) 

where the (+) sign refers to a bosonic external line and the (— ) sign to a fermionic one. 
This equation can be obtained by reversing the sign of all momenta (external as well as 
internal) and using 

A Q (-P) = A a (P) (24) 
l+Vpn M (-Po) = -(^+Vpn [7lp] (po)). (25) 

Eqs. (21) and (23) can be shown to hold at any loop order. Let us remark as a consequence 
that, for massless fields, eq. (23) holds true for the full propagators i.e., for propagators 
after the geometrical series of self energy corrections has been summed. 
The other two point functions satisfy the obvious relation 

r a a(Q) = o (26) 

which expresses the fact that a particle propagating forward (backward) in time cannot 
turn into a particle propagating backward (forward) in time by self-interactions. Techni- 
cally, eq. (26) is satisfied because its integrand is a sum of terms having poles only one 
side of the real axis with no statistical weights attached to them. 
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Turning now to the 3-point function with the momenta as indicated in fig. 4 (P + Q + 
R = 0) we can derive by the same method as above 

r aaa (P,Q,#) =0 (27) 

which is true to all order of perturbation theory. For the other functions we find it conve- 
nient to factorize out the tree vertex expression of the 3-point function and define (recall 
that 7a/?5 contains only the scalar part of the vertex as defined in eq. (15)) 

v aP s(P,Q,R) = l aps ^py (28) 

la{36\P, Qi R) 
not all indices being identical. It can be proven 

V a , s (P,Q,R) = -e 2 J 

(l +Vin l<n]( ll0 )) (A jR (Li) - A A (Li)) A Q (L 2 ) A,(L 3 ) 
L (29) 
+ (_ + %n M(/ 20 )) (A R (L 2 ) - A A (L 2 )) Ap(L 3 ) A a (L0 

+ (i + ^n^(/ 3 o)) (A fl (L 3 ) - A A (L 3 )) A*(Lx) A^(L 2 ) . 

The symbol ft is now the appropriate one for the 3-point function. Again, the constant 
factor 1/2 can be arbitrarily changed. We observe that eq. (29) appears as a sum of tree 
amplitudes since each A# — A^ combination puts the corresponding line on shell with each 
cut line carrying a weight H + m^iko)) e(l i0 ) 5(L^-M 2 ) or r] t n^(h ) e(l i0 ) 5{Ll-M 2 ) 
if the factor 1/2 is dropped in eq. (29). A graphical representation of this is given in fig. 5. 
As for the 2-point function we derive [17,18] 

V: p5 {P,Q,R) = V m {P,Q,R) (30) 
V aP6 (-P,-Q,-R) = ±V^- 5 (P,Q,R). (31) 

Where the last equation holds true for massless fields in the QED/QCD case. The + sign 
is appropriate for a fermion-antifermion-gauge boson coupling, while the — sign is for the 
triple gauge boson coupling. 

Turning to the 4-point functions they can also be expressed in terms of tree diagrams. 
They are easily obtained from [17] with the appropriate modifications of the statistical 
factors, and taking into account the appropriate Dirac structure. 

IV. R/A functions in QED and the hard loop expansion. 

In this section we write explicitly the relevant n-point functions in QED and discuss 
Ward identities as well as the hard loop approximation. 
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After performing the loop energy integration in an n-point function, by use of the 
5-function or by a contour deformation, the temperature dependent terms take the form 



n(/ ) (A fl (L) - A A (L)) F(l , /; P, Q, . . .) 

/, .n— 3 
(27r)n _! n(u) duo dl (F(fj, I; P,Q,...) + F(-u, I; P, Q, . . .)) 



(27T) 

_ 1 /' , ' 

~ 2 



where u; = \l\ and a7 is the symbol for angular integration in n-dimensions (I = l/oo). We 
have assumed massless particles for simplicity but the following argument also holds true 
in the massive case. Keeping only relevant terms we have to evaluate in practical cases 
dimensional integrals of type 

I(p) = f dw u n ~ 3 nW(w) { — Y (33) 
J rn 

where m is a mass scale typical of the external momenta components (po,p), (qo, q)--- ( a h 
external variables are supposed to have comparable sizes). This can be re-expressed with 
the variable z = uj/T as 

rn poo i_|_ p _2e 

W-Ti-rr-* I dz— (34) 

where n = 4 — 2e. These integrals are easily expressed in terms of V and Riemann ( 
functions and we find: 

l[ p) = T 2 (-r T" e r(2 + p - 26) C(2 + p - 2e) 

ITT' . yoo ) 

&\ = (1 - 2" 1 -P+ 2e ) li p) 

Defining "soft" momenta as those with all components of 0(eT), e < 1, then a Green's 
function with soft external momenta behaves as 

/( p >~T 2 (-)P (35a) 
e 

Clearly the dominant term comes from the largest value of p, i.e., from terms in the 
integrand with the highest power of L, the loop momentum. In other words, the dominant 
behavior arises from the parts in the integrand which are leading when U ~ T, hence 
the name of hard loop approximation. This is the basis of the resummation approach of 
Braaten and Pisarski [2,3] which consists in taking account of all hard loops at any order 
of perturbation theory. Several examples are discussed below. 

Consider the massless fermion self-energy in the Feynman gauge (see fig. 6). The 
general expression is (we simplify the notation by using E a instead of S aQ ) 

/An t r i 

^ 7, (Z 5 + t) l v [( g + n B (l )) (A R (L) - A A (L)) A a (P + L) 



+ {\- ™ F (po + lo)) (A R (P + L)- A a (P + L)) A S (L) 



(36) 
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For a soft fermion (p ~ eT) , we can, according to the above discussion, neglect f 3 compared 
to t> and with a change of variable in the second term (P + L — > —L) we arrive at 



-zE a (P) = -2(l-e)e 5 



2 f d n ~ Y L n B (u) + n F (u) 



(27T) 



n-1 



2w 



[£ A a (P + L)+£' A a (P + L')] (37) 



with L = (a;, /), V = (—a;, /). Consider now A Q (P + L) for example, 



A a (P + L) 



2|p + /t 



p° + oj - \p + l\ + ie a p° + U + \p + l\ + ie a 



(38) 



(we use the convention sr = e and sa = —£, £ > 0). Under the assumption that P is soft, 
we can safely neglect the second term, which behaves as 1/cu 2 , and keep only the first one 
which for large u> reduces to 



A a (P + L) 



1 



2uJ p -pl + ie a 



I 



(39) 



In the hard loop approximation, only one of the two poles of the R/A propagator is 
relevant, namely the pole associated to Landau damping [3]. The pole associated to particle 
production which is the only one at temperature, is suppressed by powers of the coupling 
constant (or equivalently by factors in 1/T). The fermion self-energy then takes the simple 
form 



t 



PL + ie c 



(40) 



where we have introduced the light-like vector L = (1,1). The dimensional and the angular 
part of the loop integration factorize, as is well-known [3,4] and we can write our final result: 



s »< p >=« £ »y 2(2^ 



2e 



(41) 



PL + ie a 

In the 4-dimensional limit, the mass term m^ h (e) reduces to the thermal mass squared [22] 

e 2 T 2 



™ 2 th = 
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(42) 



a well-known result. Although not necessary here, we have worked in n-dimensions for 
later purposes. We immediately have the useful property 



E Q (-P) = -E a (P). 
9 



(43) 



The effective propagator is defined as 



= . QP-S Q (P)) 1 } 

1 Dl(P) 

where, in the second form, -D^(P) is a scalar which need not be specified at this point. 
If P is soft, of O(eT), eq. (41) immediately shows that £ a (P) is also of C(eT), hence 
the necessity of using the resummed propagator eq. (44) rather than the bare one for a 
consistent calculation. If on the other hand P is hard, of 0(T), it appears that E a (P) 
is of 0(e 2 T) from eq. (41) and therefore the self energy gives a (negligible) correction of 
0(e 2 ) to the bare propagator. In general, hard propagators need not be resummed in a 
leading order calculation [3] . We recall the antisymmetry property of the effective fermion 
propagator which is a consequence of eq. (43) 

*S a (-P) = -*S„(P). (44a) 

It will be used later to simplify calculations. 

We could calculate, using the same technique, the vacuum polarization diagram and 
obtain the scalar 7r^(Q), 7t^(Q), i.e., the transverse and longitudinal polarization functions 
[23] which are the analytic continuations 7r L ((/ + is a , q), 7r T (qo + ie a , q) of the imaginary 
time approach. Since we do not need these results in the following we do not go into details 
and turn now to the three-point function. 

For the sake of completeness, we consider temporarily the case of massive fermions 
and work in a general gauge where the photon propagators is denoted P up {Li) A Q (Li). 
We also consider the complete expression not assuming the hard loop approximation. It 
is easy to prove (in fig. 4 we assume Q to be the photon momentum and P the incoming 
fermion momentum): 

/jn j 
7, ih + M) 7 „ (U + M) 7p ^(Li) 



(45) 



{\ +n B (h )) e(l w ) 5{L\) A a (L 2 ) A~ S (L 3 ) 
+ {\- n F (l 20 )) e(l 20 ) S(L 2 2 - M 2 ) A P (L 3 ) A 5 (Li) 
+ {\~ n F (ho)) e(ho) S(L 2 - M 2 ) Aj(Li) A (L 2 ) 

In order to prove that the usual Ward identity holds true at finite temperature we construct 
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(46) 



the scalar Q^V^^P, Q, R). Making use of the following identities: 

{U + M)@ (U + M) 5(L 2 - M 2 ) A (L 3 ) = i{U + M) 5(L 2 - M 2 ) 
(U + M)Q (U + M) 5(L 2 - M 2 ) A /5 (L 2 ) = + M) 5(L 2 - M 2 ) 

(U + M) <fy (U + M) 8{L\) A Q (L 2 ) A- 5 (L 3 ) 

= i {(U + M) A Q (L 2 ) - (^ 3 + M) A- 5 (L 3 )) 5(Lj) 

we easily derive the following identity 

QvVfaiP, Q, R) = E a (P) - X- s (-R). (47) 

This holds true not only for the e 2 T 2 terms but also for the e 2 T and of course the constant 
(T = 0) pieces. In the derivation, the thermal factors do not play any particular role. 
Knowing that, at T = 0, eq. (47) is satisfied to any loop order it can presumably be 
proven, by similar methods, that it is also true at finite temperature [24]. 

Let us now turn to the hard loop approximation. As is the case of the self-energy it 
is found that only the Landau damping contribution has to be considered. We are then 
justified in writing the propagators 

A a (L + P)= 1 



A- 5 (L-R) 



2- PL + ie a 
2w RL + ie s 7 



and recombining terms proportional to (| - np) in eq. (46) we have (massless theory, 
Feynman gauge) 

1 C rl n T 

V» p5 (P,Q,R) = --(l-e)e 2 y j^y^l 7 M t (n B (l ) + n F (l )) e(l ) 5(L 2 ) 

1 1 

(49) 



PL + ie a RL + ie$ 



i / \ f 1 dL 

~ m ^ e) J 2j^y 



dL D*t 

2e {PL + ie a )(RL + ie s ) 



The interesting feature about this formula is its independence on the retarded or advanced 
prescription on the photon momentum Q: the analytic structure of the vertex is entirely 
given by the prescription on the fermion momenta. As a consequence, in the hard loop 
approximation we can write 

V&uiP, Q, R) - V: A5 (P, Q,R) = (50) 

(not all indices being equal in the above expression). Mathematically this relation ex- 
presses the fact that the function V^ S (P,Q,R) of three variables P,Q,R, considered to 
be independent has no singularity when crossing the real qo axis. In other words 

DIscq V^ RS (P,Q,R) = (50a) 
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The independence of the variables P, Q, R is understood in the following sense: the retarded 
or advanced prescriptions on the momenta P and R respectively are not dependent on 
that of Q since (3 takes the values R, A without changing the indices a and 6. In the 
diagrammatic decomposition of fig. 5, eq. (49) is entirely represented by term a) with 
only the internal photon line being cut and the is prescription entirely carried by the 
internal fermion lines. This property results from a rearrangement of terms in the integrand 
and, as a consequence, the statistical weight attached to the cut photon line is not only 
27r(i + n B (l )) e(/ ) S(L 2 ) as in the general case but rather 27r(n B (/ ) +n F (l )) e(/ ) S(L 2 ) 
for the hard loop case. Since the angular dependence factorizes out in the integrand we can 
perform the integration over the energy and the length of the internal momentum variable 
and thus recover the thermal mass factor. This property will simplify the cut-structure of 
higher loop diagrams as will be seen later. 

The contraction of the vertex function with immediately yields eq. (47), which now 
becomes a relation between Green's function evaluated in the hard loop approximation. 
Pseudo Ward identities are also obtained [4] 

P-V a p S (P,Q,R) = -X s (R) 
R-V af3S {P,Q,R) = -X a (P) 

which are only true at the leading eT level. We define now an effective vertex 

* V Q VP, Q, R)=l" + V^ S (P, Q, R). (52) 

A dimensional analysis of eq. (49) shows that if all external vertex momenta are soft, then 
the function V£p S is of 0(1) like the bare vertex and, in that case, the effective vertex 
*V£p S has to be used for a consistent calculation. If, on the contrary, one (and therefore at 
least two) external momenta are hard then the loop correction is down by at least a factor 
e 2 T compared to the tree vertex. 



We turn now to the 4-point function and, as an example, we consider the case with 2 
external photons and an fermion-antifermion pair. We restrict ourselves, for the moment, 
to functions of type C^ RRA (P, Qi, Q2, R) and cyclic permutations on the R/A indices. 
They are the Fourier transforms of the retarded products of fields [14]. Defining the sum 
of diagrams of fig. 7 as ie 2 C£p^ s (P, Qi, Q2, R) we derive the hard loop expression (the 
superscripts \i and v refer to the Dirac indices), 

1 dL £ „ * 2 1 1 



C&siP, Qu Q2, R) = m 2 th (e) J lj^r e V % 



PL + ie a RL + is* 

(53) 

1 1 ' 

+ 



(P + Qx)L + z(e a + ep) (P + Q 2 )L + i(e a + e s ) 

when the sign of the ie terms is entirely determined by the prescriptions on the external 
legs. The Ward identities are easily obtained (we recall that one index is equal to A and 
all the others are equal to R) 

Qi^Z l5 {P, Q u Q 2 , R) = Vr. ia+P) (R, Q 2 , -R - Q 2 ) - C(/W P ' Q* ~ P ~ Q*) ( 54 ) 
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where the symbol ■ in the indices means that the corresponding R/A index need not be 
specified. 

We turn now to the case of C AAAR (P,Qi,Q 2 ,R) where we consider also the cyclic 
permutations on the indices. For this purpose it is useful to introduce the normalized 
functions C^ AAR (P, Qi, Q 2 , R) defined by [18], [20] 

C vu f PQ q R) Jn F ( P o)) s -Hn B (qio)) S ^(n B (q2o)) S - A (n F (ro)) SsA 
a(3lS ' U 2 ' n(S a APo + S ( 3Aqw + S 1 Aq20 + SsAro) (55) 

which in the hard loop approximation is given by eq. (53). This can be seen either by 
direct calculation or by using the general relation [18] (only one index a, j3, 7 or 5 is equal 
to R) 

C^P, Qu Q2, R) = {C^- S (P, Qu Q2, R)T (56) 

The function S (P, Qi, Q2, R) therefore satisfies eq. (53) which holds true then when- 
ever three of the indices are identical. We can go back to the full Green functions and we 
get, for example, 

Q^^rrra(^ Qi> Qii R) = ^arr(Ri Q2, -R - Q2) ~ T RRA (P, Q 2 , -P - Q2) 

Qi^C AAAR (P, Qi, Q 2 , R) = - — (Po) (?io) r RAA (R, Q 2 , -R - Q 2 ) /r 7 \ 

n [Po + qio) \°') 

, n B (q 10 ) n F (p + q 20 ) , , 

n {po + qio + q2o) 

We note that in the second case there appears a pre-factor in front of the vertex function 
depending explicitely on the photon momentum Qi: this is necessary since the 4-point 
function carries such factors whereas in the 3-point functions momentum Qi does not 
appear. This suggests that, in general, Ward identities take a simpler form in terms of 
the normalized functions of type and C£p g , which are analytic continuations of the 
imaginary time formalism, than in terms of the full R/A functions. This can be verified 
when considering the more complicated case of C ARAR (P,Qi,Q 2 ,R) which in the hard 
loop approximation reduces to 

CZ tAR (P,Q 1 ,Q2,R) = m 2 th (e) J ^j^^ & L v % 1 1 

1 - n F (p ) - n F (r + q 20 ) 1 + n B (q 20 ) - n F (-r - q 20 ) 1 - n F (p ) + n B (q 2 ) 

(P + Q^L + ie (P + QJL-ie ' (P + Q 2 )L-ie 

(58) 

It leads to the following Ward identity: 

Qiv-Carar(Pi Qi) $2) R) = F AAR (P, Q2, -P - Q2) - F RAA (R,Q 2 , -R - Q 2 ) 

n F (p ) n F (r + q 20 ) (59) 

H jr, x J- raa\^ V2, -a - H2) 

n {-qio) 
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PL -is RL + ie 



The last two terms appear because the prescription on the energy variable ro + (/20 is n °t 
defined by the external conditions and in such a case the R/A formalism selects a particular 
linear combination of the retarded and advanced continuations. 



V. 2-point function in the multi-loop approximation 



In the framework of the Braaten-Pisarski resummation we need the expression of 
Green's functions beyond the one-loop approximation. In particular, one is led to evaluate 
two-point functions with vertex corrections such as shown in Fig. 8a. We derive a general 
expression for the two-point function at the multi-loop level, independently of the hard 
loop approximation to which we return at the end of the section. 

In the case of interactions involving only three particle vertices, it is always possible 
to cut the diagram in such a way as to have only two particle intermediate states (Fig. 8b). 
This allows us to express the self-energy in terms of the vertex functions Y a ps 



(60) 



where the various Y functions carry the superscripts (j), (k), and (j + k+ 1) to denote 
the number of loops at which they have been respectively evaluated. Introducing the 
normalized functions as in eq. (28) and regrouping terms we obtain (dropping irrelevant 
spin factors) 

W V«> J ( 27r )n 

i + W p] (Po)) Di scp [A fl (P) V$ A (P,Q,-R)A R (R) vj» R (-P, -Q, R)] (61) 
+ Q + W* ] (ro)) Dis^ [A A (P) V$ a (P, Q, -R)A r (R)V$ r (-P, -Q, R)] 

where the symbol Discp applied on a function F a ps(P, Q, R) means taking the discontinuity 
in the energy variable po of the function as defined by 

Discp F R/35 {P, Q, R) = F Rf3S {P, Q, R) - F A/3S (P, Q, R) (62) 

(see the discussion around eq. (50)). The ingredients to derive this relation are the same as 
those to obtain the one loop expression, namely the diagonality of propagators, the validity 
of eq. (27) for the dressed vertex and the property that terms with poles only on one side of 
the real axis of the loop energy variables give a vanishing contribution when no statistical 
weight, depending on the loop energy variable, is attached to them. Despite a rather 
cumbersome notation, the structure of eq. (60) is rather simple and similar to the one-loop 
result. Consider, as an example, the calculation of Y RR (Q). The function in the first 
line, whose Discp should be evaluated, is trivially obtained with propagators and vertex 
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functions in the retarded P and Q momenta (recall also the crossing relation eq. (13)). 
Momentum R is necessarily retarded because of momentum conservation. Likewise the 
second line is constructed from the retarded momenta Q and R. However, P is now of the 
advanced type since, using po = ro — qo and recalling that taking the discontinuity in ro 
puts ro as the real axis, the prescription for po is that of —qo, i.e., advanced. Taking the 
discontinuity of these products of propagators and vertices generates two types of terms: 
Discp A R (P) = 2ir e{p ) 5(P 2 - M 2 ) puts the P line on shell while Disc P V R/3S {P, Q, -R) 
amounts to taking the cut contribution of the vertex. 
Likewise it is easy to prove 

rflg(Q) = o. 

As an application we may consider the case when V a /3$ consists in a one loop diagram 
and study Discp V a p$(P, Q, R). We have already seen that the vertex is expressed as a sum 
of the diagrams (eq. (29) and fig. 5). Taking the discontinuity with respect to po means, 
according to eq. (62), putting the internal momentum lines carrying the index a on shell. 
An examination of eq. (29), leads then to taking the line L 2 = Li + P on shell when L\ 
is cut and L\ = L 2 — P on shell when L 2 is cut. This we symbolize by a cross on the 
corresponding line (see fig. 9). We use a different symbol to denote this cut because the 
statistical weight associated to taking Discp is (\ + i]pn^ np ^(po)) rather that the thermal 
factor appropriate for the internal line. Therefore, Discp V a ps in eq. (61) picks up the 
2-particle cut contribution in p while Discp A R (P) selects the pole in p Q . One can easily 
convince one-self, by such arguments that the self-energy in the multi-loop approximation 
retains a tree structure as it does at one loop. 

In the next section we will be led to consider the imaginary part of the two point 
functions which is evaluated by constructing T RR (Q) — T AA (Q). It can be written as 

r RR (Q)-T AA (Q) = -te 2 J ^| Q+W^kpo)) 

DiscpAp(P)VWP, Q, -R)A R (R)V AAR (-P, -Q, R) 

- Disc P A R (P)V RAR (P, Q, -R)A A (R)V ARA (-P, -Q, R) 
+ (l + VRn [m] (r )) 

Discp A A (P)V ARA (P, Q, -R)A R (R)V RAR (-P, -Q, R) 
-DiscpAp(P) VW((P, Q, -R)A R (R)V ARR (-P, -Q, R) 

The difference of the Discp or Discp expressions above have a simple interpretation. Let 
us remark that in the first square brackets the two terms differ by the indices associated 
to momenta Q and R while in the second ones the terms differ by the indices associated 
to momenta Q and P. Consider the Discp case. As we have seen, taking the discontinuity 
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with respect to the variable po amount to putting the energy p on the real axis. By 
"momentum conservation", tq = —po — qo, the retarded/ advanced prescription on the 
internal momentum R is now entirely determined by that on the external momentum Q: 
ro, in that sense, becomes a function of qo. This is expressed by the fact that the Discp 
terms differ by their indices in both Q and R. Keeping this in mind, the square brackets 
then isolate the discontinuity in qo of the considered functions and we may write 

DiscQ Discp A R (P) V RRA (P,Q,-R) A R (R) V AAR (-P,-Q,R) 

= Discp A R (P) V RRA (P,Q,-R) A R (R)V AAR (-P,-Q,R) (64) 

-Discp A R (P) V RAR (P,Q,-R) A A (R)V ARA (-P,-Q,R) 

where R is considered as a function of Q, when Discp is evaluated. 

If we turn now to the hard loop approximation we find that the singularity struc- 
ture of the effective vertices simplifies considerably. A case in point is the QED vertex 
studied in the previous section where it was shown that (in the hard loop approximation) 
DiscQ V a Rs{P, Q, R) = (eq. (50a)). This function enters the calculation of the vacuum 
polarization tensor TTp U (Q) to be considered shortly. In eq. (64), as a consequence, we 
can write Discp Discp for the first term and Discp Discp for the second one. Further- 
more, from eq. (49) it appears that the dependence on P and R factorizes in the integrand 

which we can write as a product of two functions f a (P) f$(R): the double discontinuity 
becomes then a product of discontinuities so that 

r RR (Q) -T AA (Q) = -te 2 j ^JL^l-n F (p Q )) Discp/jftp) Discp/f (P) 

+ ((1 - n F (r )) Discp/i 1} (P) Discp/^P) 

(65) 

where we have considered the case with an internal fermion loop. Combining both terms 
and using the detailed balance relation 

n F (r ) - n F (p ) = (1 - e^°)n F (r )n F \-p ) 

it comes out 

/Jn— 1 p r 
T^j— r (1 - e^) J dp dro n F (r )n F (-p ) S(q + p - r ) 

^Discp fS\p) Discp f R 2) (R). 

(66) 

This is to be compared to the formula in [5], [6] which expresses the discontinuity of the 
two-point function as an integral over real energies 

Jn-l 



^I^Disc T / (1) (PO, P)f (2) (d° -P°,Q-P) 

^ ' Tin 



Po 

I 



J dwdw'n F (u')n F (uj')5( qo - U - u/)2irp 1 (u; 1 p)p 2 (u/ 1 q - p) 

(67) 
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where the spectral density pi is defined as, for example, pi(u,p) = 27rDisc f^ l \po,p)- 

The structure of eqs. (61), (63) is not changed, if instead of the bare propagators 
A a (P), As(R) we use in these equations the effective propagators *A a (P), *Ag(R). The 
poles near the real axis are shifted away from it and become singularities in the same half- 
plane of the loop energy variable. In particular, the crucial ingredient in deriving these 
equations, namely that an advanced propagator has no pole in the lower half-plane and 
a retarded one has no pole in the upper half-plane, is not affected. When evaluating the 
discontinuities of these effective propagators we will not only get the pole contributions as 
in the bare case but also two particle cuts associated to the Landau damping mechanism. 



VI. Real soft photon production in a quark-gluon plasma. 



The problem of photon emission in a quark-gluon plasma has already been considered 
for several cases: soft virtual photon at rest [5] or moving [6], [25] and hard real photon 
[7] , [8] . We apply the above formalism to soft real photon [26] . We assume massless quarks 
and we introduce the strong interactions coupling, denoted g, assuming j< 1. The photon 
production rate [27], assuming qo positive for simplicity, 

Q °^Fq = ~ (27r) 3 ?lB( ^ IlrLl1 ^ (^) I Retarded (68) 

is related to the trace of the polarization tensor which in the following we denote for short 
Rr(Q), after summing over the photon polarization states. 

Consider first the one-loop contribution to Im H R (Q) before resummation, i.e. using 
bare propagators and vertices. Using eq. (18), it is trivially found that Im H R (Q) vanishes 
(but Re H R (Q) 7^ 0). The reason for this is simple: the only kinematical configurations 
possibly contributing to Im IT^ are the collinear decay of the photon into a qq pair or the 
collinear emission or absorbtion of the photon by a quark or an antiquark in the plasma. 
However helicity conservation at the 'yqq vertex forbids these processes. 

We turn now to the effective theory and evaluate the same "one loop" diagram using 
effective propagators and vertices as shown in Fig. 10 a). The tadpole diagram (Fig. 10 
b) vanishes because it is traceless [5]. In the framework of the resummed perturbative 
series, the self-energy like diagram gives the dominant contribution to the production of 
soft virtual photons. We have 



/Jn p 



DiscpTr 



*s R (pyv£ RA (p, q, -Rys R (Ryvz RA (p, Q, -R) 



(69) 



Some words of explanation are required concerning this equation. The effective prop- 
agator *S a and vertex *V ( ^ 3(5 are understood to contain QCD corrections and not QED 
corrections. At the order at which we do the calculation it simply amounts to substituting 
e 2 — > CVg 2 in eq. (42) to take into account the change in coupling as well as the color 
factor. We thus define the thermal mass to be now 

m 2 th = C F ^—, (70) 
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all other equations in the previous section being unchanged. We have applied the property 
of eq. (31) to the effective vertex eq. (52) so that the same vertex function appears twice 
in the integrand above. Finally, the expected term proportional to (| — ni?(ro)) can be 
reduced to the one above after the change of variable ro = — po and the use of eqs. (31) 
and (43), hence the factor 1 — 2 n F (po). Writing out the explicit form of *V^ /3S we thus 
have to calculate the diagrams of Fig. 11. We ignore the first one which does not present 
any difficulty and turn to the second (or equivalently the third) one which will be shown 
to present a collinear divergence. We have: 



-*n*(Q)|6 = -eV th {e) J^J Ij^-Al - 2 n F (p )) 

Tr(*S R (P) I *S R (R) I) 



(71) 



Discp- 



(PL + ie)(RL + ie) 



We choose to carry the f dpo integration by closing the contour in the upper half-plane. 
Writing Discp explicitely as the difference of two terms, 

Digc Tr = Tr(*S R (P) I *S R (R) I) _ Tr(*S A (P) j *S R (R) j) 
P (PL + ie)(RL + ie) (PL - ie)(RL + ie) 

we see that no contribution arises from the first term since all its singularities are located 
below the real axis. On the contrary, the second one contains a pole in the upper half plane 
(1/(PL — ze), coming from the hard loop in the effective vertex as well as a singularity 
from the resummed propagator *Sa(P)- Let us concentrate on the pole contribution. It 
comes out 

/d n P 1 C df 
(2^2 j WF^ 6(PL) {1 ~ 2nF{Po)) 
Tr(*S A (P) I *S R (R) I) (73) 
(QL + ie) 

where the po integration is understood now to run along the real axis. The denominator 
QL + ie is nothing but RL + ie where the 5-function constraint has been used. Both Q and 
L being light-like this factor leads to a collinear divergence when the angular integration 
j dL is performed. From now on, we are only interested in this diverging part neglecting 
all finite terms in the calculation. To evaluate the residue at the pole, it is enough to set 
L = Q = Q/q in the above, except of course, in the denominator, leading to the rather 
simple expression 

/d n P * 
S(PQ) (1 - 2n F (p ))TtCS A (P) £ *S R (R) #) 



i r di 

2j 



dL 1 

" 2e QL + ie 

(74) 



18 



The angular integration is understood in n — 2 dimensions (n = 4 — 2e) and its real 
part has a pole in e. Therefore we define 



7(e) 1 [ dL q 



e 2 J (2ny~^QL + ie 

neglecting a finite imaginary piece. The diverging contribution from the two diagrams 
containing one hard loop effective vertex is therefore 



n R (Q)\ b+ csin 9 = 2 e 2 HhMM I _*^ S (pQ ){1 _ 2n F ( Po )) 

Tr(*S A (P) ® *S R (R) ®) 



(76) 



To sum up, the diverging contribution arises from the momentum configuration where 
both PL and RL vanish that is when the fermion propagators coupling to the external 
photon both approach their mass-shell condition. It is interesting and somewhat paradox- 
ical that such a divergence is a property of the resummed theory which tells us to use the 
effective propagators *A a (P) and *A$(R) rather than the bare ones. Had we used the 
latter, coming back to eq. (71), we would have found that the trace reduces to 

Trflftl = 8PL RL (77) 

cancelling both poles at the origin of the collinear divergence. 

It is worth noting that the collinear divergence in eq. (75) is related to the vanishing 
mass of the photon. If Q 2 were slightly off-shell the pole in eq. (75) would never be reached 
in the time like case or would be defined through a principal value prescription in the space 
like case. 

Let us turn now to the last diagram with two hard loop effective vertices. It is 
explicitely 



-m R (Q)\ d = -fmi h (e) J^IJ \ J 



dL, 



Discp ( 



(27T) 1 - 26 

L X L 2 Tr(*S R (P) ti *S R (R) l 2 ) 
V (PL 1 + ie) (R^ + ie) (PL 2 + ie) (RL 2 + ie) ■ 



(78) 



As in the previous case, a collinear divergence will arise from the collinear configuration 
L\ = Q, leading to a pole in e. It is to be noted that the configuration L 2 = Q, together 
with the collinearity condition on L\ does not lead to a double pole because of the L\L 2 
factor. Setting thus L\ = Q in the denominator we construct in fact a combination of 
type Q nV^g(Pj Qi —R) which can be immediately reduced via the Ward identity eq. (47). 
More precisely we have in the integrand 



™ 2 th 



( C )I / - dL * — , 9Llh, = I(S Q (P) - E 5 (R)). (79) 

y, 2j (27T) 1 "* (PL 2 + ie a )(RL 2 + ze s ) <Z 1 ) 1 " 1 ) 
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Recasting the f dpo integration as in eq. (73) and exhibiting the collinear divergence 
it comes out 

-in R \ d , sing = -2ie 2 m 2 th (e) ^ f -f^ 5(PQ) (1 - 2 n F (p Q )) 

C J (27r) (80) 
- Tr(*S A (P) fa *S R (R) (E A (P) - E R (R))) 

The factor 2 is introduced to account for the case when L 2 is collinear to Q. Using 
the definition eq. (44) of the effective propagator the trace term is reduced to 

% (Tr(* S A (P)fa) ~ Tr(*S R (R)fa)) - q Tr(* S A (P) fa *S R (R) fa) (81) 

It is immediately apparent that the last term exactly compensates the single effective 
vertex pole contribution eq. (76). We are thus left, adding all pieces together, with 

n s (Q)U 9 = i 2 e ' 2fe> | ^r*(«5)(i - ^( Po) ) (g2) 

{Tr(*S A (P) ff)-Tr('S R (R) ($)) 

The Ua(Q) 2-point function is easily deduced from this equation by replacing *S a (P) 
by *Sa(P) and the same for *S$(R). We introduce now the usual parametrization of the 
effective fermion propagator 



(83) 



where P s is the light-like vector P s = (l,sp) and D S (P) is a scalar function whose ex- 
pression can be found in [5], [6]. The sum runs over the two propagating modes of the 
thermalized fermion. From the general property of 2-point functions, eq. (21), and the 
definition eq. (44), we can write 

*S A (P) = -fg fi (f)r = iV m -. Ps — (84) 

1 ' y y " ^ (D R (p + ie,p))* K ' 

which expresses the advanced effective propagator as the complex-conjugate of the retarded 
one. It is customary to use the notation 

= a a {P)^mP 8 {P) (85) 



R,A 



We can now construct the imaginary part of the polarization function 
ImU R (Q) = 2 e 2 ^ M J S(PQ) (1 - 2n F (p )) 

it (ps(P)Trf s fa + p s (R)Trft s 

s 
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(86) 



The traces are easily evaluated and give 4(1 — spq) and 4(1 — srq) respectively, which 
when taking into account the 5-function constraint in the integrand reduced to 4(1 — spo/p) 
and 4(1 — stq/t) leading to 

ImlMQ) = 8 e 2 ^MM J S(PQ) (1 - 2n"(po)) 

^((i-^) A (P) + (i-£r5) ft(B )) (87) 

s ^ 

If we remember that the j3 s function are proportional to the same (1 — spo/p) type 
factor with a positive coefficient we see that the integral does not vanish and therefore 
the soft fermion loop contribution to the 7 production rate in a plasma is divergent. It is 
interesting to note that the condition 5(PQ) enforces that P (and also R since S(RQ) = 
5(PQ)) be space like (\po\ < p), (\fo\ < r). The residue of the collinear pole is entirely 
due to the Landau damping contribution to the fermion effective propagator. 

An alternative derivation of eq. (87) has recently been given in the imaginary time 
formalism [26]. It seems appropriate now to contrast the two approaches. Since we work 
here in the real time formalism we can use the usual Dirac algebra familiar from the T = 
case and we do not have the added complication to redefine the algebra in Euclidean 
space. A more important point derives from the fact that we calculate the full retarded (or 
advanced) 2-point function and not only its imaginary part. We manipulate 4-dimensional, 
or rather n— dimensional integrals and the integrands, except for statistical weights, keep 
a covariant form very similar to the T = case (see e.g., eqs. (71) and (78)): because we 
work in n— dimensions we can evaluate the integral over the energy by closing contours 
in the appropriate half-planes, as we do in eq. (71) and eq. (78), to pick-up only the 
pole contributions which lead to the singular behavior; since after this manipulation the 
integrand keeps its covariant form we easily use the Ward identity (eq. (47)) to extract the 
final result. In contrast, in the ITF approach, the imaginary part of the 2-point function 
is directly calculated as an integral of spectral functions over the real energy axis: the 
various pieces of the integrand have to be decomposed into their principal values and 
5— functions to extract the divergent behavior and since the covariant form is lost in favor 
of the spectral functions the Ward identity cannot be directly used to show the partial 
compensation between the different terms. An explicit calculation of all terms is necessary 
to obtain the final results. 

VII. Conclusions 

In this work, we have pursued the study of the R/A formalism and derived, in par- 
ticular, an important equation expressing the discontinuity of the 2-point function at any 
loop order. We have also formulated, in this approach, the hard loop expansion for the 
case of QED. We used the R/A formalism to calculate the rate of real soft phton produc- 
tion in a QCD plasma and found the result (eq. (87)) to be divergent in agreement with 
a recent calculation performed in the imaginary time formalism [26]. Compared to that 
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work we find the calculation in the R/A formalism somewhat simpler, since, using contour 
integration, we can extract the relevant contributions of poles in the complex energy plane. 
We avoid distinguishing between principal part and 5-function contributions of a pole as 
when the integration is performed on the real axis. Also the cancellation of the single hard 
vertex loop diagrams with the double hard vertex loop diagram is clearly attributed to the 
thermal Ward identity eq. (47). The structure of the remaining diverging term involves 
the spectral function of the effective soft fermion propagator. It is interesting to see how, 
in the framework of the effective theory diagrams with different number of loops partially 
compensate one another. 

There are several ways the divergent result of eq. (87) could be regularized. In [26] 
it was proposed to introduce a soft cut-off of order gT. Another possibility is to include 
corrections to the hard internal fermion lines which are the cause of the collinear divergence 
when they approach the mass-shell condition. This is in the spirit of several attempts to 
define the damping rate of a hard fermion. For example, including a width of 0(g 2 T) 
on the hard internal fermion lines [28] would displace the poles away from the real axis 
and replace the 1/e divergence by ln(l/gr). Finally, we have considered in our calculation 
only the contribution of diagrams with two internal soft fermion lines (see Fig. 11). As 
discussed above the hard fermion loop is kinematically allowed but it does not contribute 
at the lowest order using bare propagators because of helicity conservation at the vertices. 
Introducing corrections to hard loops avoids this constraint. It would be interesting to 
compare the order of such contributions to the result of eq. (87) . 
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Figure Captions 



Fig. 1 The general "real-time" contour in the complex time plane. 

Fig. 2 The QED vertex in the R/A formalism : a) with an incoming fermion and an incoming 

anti-fermion; b) with an incoming fermion and an outgoing fermion. 
Fig. 3 The generic two-point function. 
Fig. 4 The generic three-point function. 

Fig. 5 The three-point function as a sum of tree diagrams. The / on a line means that the 

line is put on mass-shell and carries a statistical factor: e(ko)(^ + rjin^ Vi \lio))S(Lf). 
Fig. 6 The fermion self-energy diagram. 

Fig. 7 The photon-photon-quark-antiquark four-point function in QED. 

Fig. 8 A diagram in the effective theory: the two-point function with vertex corrections; a) 

one loop vertices; b) dressed vertices. 
Fig. 9 Taking the discontinuity with respect to po of a vertex up to one loop. The cross 

represents the action of taking the discontinuity (i.e. putting the corresponding line 

on shell). 

Fig. 10 The soft fermion loop contribution to the rate of soft real photon production in a 
QCD plasma. 

Fig. 11 The same as Fig. 10 displaying the structure of the effective vertices. 
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